CALT-68-1783 
UCSD/PTH 92-17 
DOE RESEARCH AND 
DEVELOPMENT REPORT 



Baryons Containing a Heavy Quark as Solitons 
(REVISED VERSION) 

Elizabeth Jenkins and Aneesh V. Manohar 

Department of Physics, University of California, San Diego, 
9500 Oilman Drive, La Jolla, CA 92093-0319 



Mark B. Wise 

California Institute of Technology, Pasadena, CA 91125 

Abstract 

The possibility of interpreting baryons containing a single heavy quark as bound 
states of solitons (that arise in the nonlinear sigma model) and heavy mesons is 
explored. Particular attention is paid to the parity of the bound states and to the 
role of heavy quark symmetry. 
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Many properties of baryons B containing light u and d quarks suggest that they 
can be represented as sohtons. For example, they have a mass of order 1/a and 
the cross section for e+e" — > BB is of order e"-*^/", where a — 1/Nc is taken as a 
small quantity.^ Skyrme originally suggested that baryons are solitons in the chiral 
Lagrangian (i.e. nonlinear sigma model) used to describe pion self interactions.^ 
The solitons of the chiral Lagrangian have the right quantum numbers to be QCD 
baryons,^ provided one includes the Wess-Zumino term."^ The model of QCD baryons 
as solitons has been used to compute many of their properties. ^'^ The same large Nc 
power counting that suggests that baryons containing only light quarks are solitons, 
also suggests that baryons containing a single heavy quark Q {mq ^ ^qcd) can be 
represented as solitons. 

Recently there has been considerable progress in understanding the properties of 
hadrons containing a single heavy quark. In the limit uiq oo, QCD has a heavy 
quark spin-flavor symmetry^'® that determines many properties of hadrons containing 
a single heavy quark. In the mq — > oo limit, the total angular momentum of the light 
degrees of freedom (i.e. hght quarks and gluons), 

Si^S-Sq, (1) 

is conserved. Consequently mesons and baryons containing a single heavy quark can 
be labeled by s^, and provided 7^ 0, they come in degenerate doublets^ with total 
spins 

s± = s^±l/2 (2) 

formed by combining the spin of the heavy quark with the angular momentum of the 
light degrees of freedom. 

Callan and Klebanov originally suggested an interpretation of baryons containing 
a heavy quark as bound states of solitons of the pion chiral Lagrangian with mesons 
containing a heavy quark. ■'^^'■^■'^ In this paper we examine this model for heavy baryons. 
An important aspect of this work is that the consequences of heavy quark symmetry 
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are taken into account. Heavy quark symmetry implies that a doublet of mesons con- 
taining the heavy quark must be considered and that both members of the doublet 
can occur in the bound state. In this regard the work presented here differs from that 
of Callan and Klebanov. An important feature of the work of Callan and Klebanov 
was that the ground state of the heavy baryon is in an L = 1 partial wave, and thus 
the lowest mass baryon containing a heavy quark has positive parity in accord with 
experiment. In this paper both the soliton and meson are taken as infinitely heavy 
and it is assumed that the soliton-meson potential is minimized at the origin. Conse- 
quently, the spatial wavefunction for the ground state is a Dirac delta function S^{x) 
and corresponds to the lowest energy classical configuration which has the soliton and 
heavy meson located at the same spatial point. We give a simple explanation of the 
Callan-Klebanov parity fiip appropriate to this description of the bound state. 

The lowest lying heavy mesons with Qqa (qi — u, q2 — d) fiavor quantum numbers 
have S£ — 1/2 and come in a doublet containing spin-zero and spin-one mesons. For 
Q — c these are the D and D* mesons. The interactions of these heavy mesons with 
pions are described by a chiral Lagrangian which is invariant under both heavy quark 
spin symmetry and chiral SU{2)l x SU{2)ji symmetry* Chiral SU{2)l x SU{2)ji 
symmetry is spontaneously broken to the vector SU{2)v subgroup, so the Gold- 
stone boson manifold is SU{2)i x SU{2)ji/SU{2)y. A general Goldstone boson field 
configuration is obtained by making a space-time dependent chiral transformation 
(giix), gR{x)) on a standard vacuum state. The Goldstone boson fields 5(x) in a 
G/H sigma model are defined by^^ 

g{x) = E{x)h{x) (3) 

where g{x) G G describes the locally rotated vacuum, h{x) e H, and E{x) is generated 
by Ug independent broken generators, where Ug is the number of Goldstone bosons. 



■k We will restrict the analysis to the case of two light flavors; the generalization to three flavors 
is straightforward. 
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Different choices for the Ug broken generators lead to different (but equivalent) non- 
linear realizations of the spontaneously broken chiral symmetry. The standard E 
basis for the QCD chiral Lagrangian is defined by 



{9l{x),9r{x)) ^ (S(a;),l) • {V^{x),V^{x)) = {E{x)V^{x),V^{x)) , 



(4) 



where the broken generators are chosen to be the SU (2) ^ generators T^. The field 
E(x) is equal to gL{x)g\^{x), and transforms under chiral SU{2)l x SU{2)ji as 



S(a;) ^ L T.{x) B} , 



(5) 



since gL{x) — > Lg^^x), gR^x) — > RgR^x). The pion fields occur in the 2x2 matrix 

-7r7V2 



M 



TT 



(6) 



where 



E = exp 



2iM\ 



(7) 



/ J ' 

and / is the pion decay constant. Experimentally / ~ 132 MeV. In the large Nc limit 
/ ~ (recall that a — l/Nc). For considering the couphngs of heavy mesons to 

the pseudo-Goldstone bosons, it is useful to introduce the ^ basis defined by 

{gL{x),gR{x)) ^ {c{^),CHx)) ■ {V^{x),V^ix)) = (cixnix) , CHx)V^{x)) (8) 

where the broken generators are chosen to be the axial generators Tf — T^. Under 
chiral SU{2)l x SU{2)ji, ^ transforms as 



^{x)^L^{x) U^x)=U{x) ^x) , 



(9) 



where [/ is a complicated function of L, R and the pion fields and in general depends 
on the space-time coordinates. However, for an unbroken SU{2)v transformation. 
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V — L — R, U{x) is equal to V, and is constant. Comparing eq. (8) with eq. (4), we 
see that ^ and E are related by 

E(x)=e'(^)- (10) 

The Goldstone boson manifold SU{2)i x SU {2)ji/SU (2)y is diffeomorphic to the 
group manifold SU{2), which is diffeomorphic to the three-sphere S"^. The standard 
identification of S with S"^ is obtained by writing 

E = a + i6-f, a2 + |6|2 = i, (11) 

and identifying S with the point (a, h) E S^. The S basis provides a single-valued 
coordinate system that covers the entire surface of SU{2), because 

{gL, QR) = (Eli:, V^) = (eV^, V^) (12) 

implies that E = E'. 

Unlike the E basis, the ^ basis is not a well-defined coordinate system on SU{2). 
Two fields ^ and ^' denote the same Goldstone boson configuration if 

{gL.gR) = (i{x)V^{x),i\x)V^{x)) = {i'{x)Vl{x),i'\x)Vl{x)) . (13) 

This equahty is satisfied if and only if 

= Hh^t heSU{2). (14) 

It is straightforward to solve this equation for h e SU{2). Denote the general SU{2) 
matrix ^ by 

^ = a + i(3-f, a2 + |/3|2 = l. (15) 

If a 7^ 0, then ^ and are equivalent. However, if a = 0, all configurations of 
the form i/3 ■ r are equivalent, irrespective of the value of /3. If the vector (a,/3) 
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is considered to be a point on the three sphere S^, then we see that points on the 
northern and southern hemisphere of obtained by inversion through the origin are 
equivalent, and the entire equator is identified to a point. (The resulting manifold is 
still equivalent to S^.) Thus the $, basis gives a coordinate system for the Goldstone 
bosons which is well-defined except near the equator of , ^ = iP ■ f. In the S basis, 
the points on the entire three sphere S'^ correspond to inequivalent field configurations, 
and the E basis gives a well-defined coordinate system on all of S^. 

The chiral Lagrangian for heavy mesons is easy to construct in the ^ basis since the 
heavy meson fields transform in a simple way under parity. Heavy quark symmetry 
is incorporated into the chiral Lagrangian by defining a heavy meson field which 
contains both pseudoscalar and vector mesons. The heavy meson field for the ground 
state Qqa mesons is written as a 4 x 4 matrix-^^ 

^a = ^^^[Pa>^-Pa75], (16) 

where is the heavy quark four- velocity, v"^ — 1. We shall work in the rest frame of 
the heavy mesons where v^^ = (1,0). The fields Pa, destroy the heavy pseudoscalar 
and vector particles that comprise the ground state S£ = 1/2 doublets. The vector 
meson field is constrained by v^P*^ = 0. Under SU (2)q heavy quark spin symmetry 

Ha ^ SHa , (17) 

where S e SU(2)q, and under chiral SU{2)l x SU{2)r 

Ha - H,UI . (18) 

It is also convenient to introduce 

Ha = I'Hh' = [P:1^^ + Phs] ^ . (19) 
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Under parity 



H{x^,x) -^Y^H{x^,-x)i 



,0 



(20) 



and 



—X 







(21) 



The chiral Lagrangian density for heavy meson-pion interactions is 



14,15,16 



(22) 



where the eUipsis denotes terms with more derivatives and repeated a, b indices are 
summed over 1,2. Eq. (22) is the most general Lagrangian density invariant under 
chiral SU (2) ^ x SU (2)^, heavy quark spin symmetry and parity. It is easy to general- 
ize this Lagrangian density to include explicit SU{2)i x SU{2)^ symmetry breaking 
from u and d quark masses and explicit SU{2)q symmetry breaking from ^qcdIthq 
effects. 

The couphng g determines the D* — > Dtt decay width, 



The present experimental hmit,^^ r(D*+ D^7r~^) ^ 72 keV, implies that g"^ ^ 0.4. 
In this work the sign of g plays an important role. Applying the Noether procedure 
to the Lagrangian density in eq. (22) and to the QCD Lagrangian density gives 



where T is a traceless 2x2 matrix and the eUipsis denotes terms containing the 
pion fields. Treating the quark fields in eq. (24) as constituent quarks and using the 



*+ 




(23) 



qaliJ,l5Tabqb = -gTrHaHb-f^-f5Tba + ■■■ , 



(24) 
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nonrelativistic quark model to evaluate the D* matrix element of the l.h.s. of eq. (24) 
gives^^ g — 1- (A similar estimate for the pion nucleon coupling gives qa — 5/3.) In 
the chiral quark model-*^^ there is a constituent-quark pion coupling. Using the chiral 
quark model gives ~ 0.75. 

The soliton solution of the SU (2) l x SU (2)ij chiral Lagrangian for pion self 
interactions is 

E = A{t)J:o{x)A-\t) , (25) 

where 

Eo = cxp {zF{\x\)x-t) , (26) 

and A{t) contains the dependence on collective coordinates associated with rotations 
and isospin transformations of the soliton solution, 

A^ao + ia-T , (27) 

where OQ+jap = 1. The soliton has the quantum numbers of a baryon containing light 
u and d quarks. For solitons with baryon number one, the function F(|x|) satisfies 
F(0) = — TT and F{oo) — 0.* Its detailed shape depends on higher derivative terms in 
the chiral Lagrangian for pion self interactions. Since F increases as \x\ goes from zero 
to infinity, it is expected that F'{0) is positive. The baryons containing u and d quarks 
have wavefunctions that are functions of the a^. For the neutron and proton states; 
\p t) = (l/7r)(ai + ia2), \p I) = - {i / n) {ao - ia^) , \n t) = {i / n) {ao + ia^) and \n |) = 
— (l/7r)(ai — ia2) ■ The soliton So(a:) is a winding number one map from spacetime 
into the Goldstone boson manifold. In the ^ basis, the soliton solution has the form 
(up to an overall sign) 

Co^exp (^F(|f|) i;-r/2). (28) 

As F varies from F = at spatial infinity to F = — tt at the origin, varies 
from ^0 — 1; the north pole of S^, to — —i^^ ■ the equator, with intermediate 

* There appears to be a sign error in the formula for the baryon number current used in Ref. 5. 
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points covering the northern hemisphere of S^. However, we have aheady seen that 
the Goldstone boson manifold in the ^ basis is the northern hemisphere of S^, with 
the equator identified to a point. Thus the sohton solution eq. (28) is a winding 
number one map. It looks singular at |x| = because the ^ coordinate system for the 
Goldstone boson manifold is singular at the equator of S^. However, this singularity 
is a reflection of the coordinate singularity on the Goldstone boson manifold in the ^ 
basis and is not a physically singular field configuration. 

The form of the Lagrangian density for heavy meson-pion interactions in eq. (22) 
is not very convenient for discussing heavy meson-soliton interactions because of the 
coordinate singularity in ^. It is convenient when dealing with matter fields interacting 
with solitons to redefine the fields so that the Lagrangian density is expressible in 
terms of S. For example, introducing new heavy meson fields and P^j* , defined by 

K = H^iba , (29) 
the SU{2)l X SU{2)ji transformation law in eq. (18) becomes 

K ^ H'^rI , (30) 
and the Lagrangian density in eq. (22) becomes 

C = -iTrH'.v^d^H', + '-TrH',Hy{T)d^T.)ba 

+ |Tr^X7V(S^5.S),, + ... (31) 

The parity transformation rule for the primed heavy meson fields is a little more 
complicated, 

H',{x\ x) ^ i^Hl{x\ -f)7° T.l{x\ -x) . (32) 

Note that in a background Goldstone boson field configuration of a sohton located 
at the same spatial point as the heavy meson the factor of becomes —1, whereas 
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= 1 for a heavy meson infinitely far from the sohton. This relative minus sign 
is the source of the parity flip that gives positive parity heavy meson-soliton bound 
states. 

In the large N^. limit, the baryon solitons B containing light u and d quarks are 
very heavy and time derivatives on the pion fields can be neglected. Consequently 
the interaction Hamiltonian 

Hi --{J Tr^^/^h'75pt9,E],, + ... , (33) 

with S given by eq. (25) determines the potential energy of a configuration with a 
baryon B at the origin and a heavy P' or P'* meson at position x. Using eq. (16) 
and taking the trace in eq. (33) gives 

Hj^gJ d'x [C^-t(f)n'(f) + i^t(f)p-^-(f) + ie^^^p-^t(^)p-fc(^)] 



X 



+ ... 

ba 



uu 

(34) 

Assuming that in attractive channels the potential energy is minimized at x — where 
the heavy meson and baryon soliton are located at the same point, it is eigenvalues 
of the potential operator at the origin that are needed to determine the spectrum of 
bound states. Using eq. (34) 

Vj{6) = gF'{0)[P;:^^Pl + P'a^K' + ^e*'^i^a"^n'*'] 

x[(ao + m ■ f)r-^(ao - m ■ r)]6o + ••• ■, (35) 

where the ellipsis denotes the contribution of terms in the chiral Lagrangian with more 
than one derivative. It is easy to diagonalize the potential (at the origin) matrix that 
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arises from taking the truncated basis of nucleon- heavy meson product states (i.e., 
\p T)|-Pi')) etc.). On this space it is straightforward to show that* 

Vi{Q) = -^^F'(0)(5/ - 3/2)(/2 - 3/2) + ... , (36) 

where Si denotes the total angular momentum vector of the light degrees of freedom 
(soliton and heavy meson combined) and / denotes the total isospin vector. The 
eigenstates of Vj{0) have definite isospin, spin and angular momentum for the light 
degrees of freedom. They are denoted by |/, s, Sf), where / is the total isospin, s the 
total spin and si is the angular momentum of the light degrees of freedom. 

The spatial wavefunctions of the eigenstates of eq. (36) are 5^ (x), and arc parity 
even. The parity of the meson-soliton bound state is also even, because the primed 
heavy meson fields are odd under parity at infinity, and so are even under parity at 
the origin, as noted below eq. (32). The unprimed heavy meson fields have a simple 
transformation law under parity, and do not have a relative minus sign between the 
parity at infinity and parity at the origin. However, in the ^ basis, the wavef unction 
of the bound state contains a factor of the form t-x near the origin. Thus the parity 
of the sohton is also even in the ^ basis because the negative parity of the spatial 
wavef unction is combined with the negative intrinsic parity of the heavy meson field. 
The factor r ■ x appears singular at the origin, but that is because of the coordinate 
singularity in the ^ basis. Any physical quantity is well-defined at the origin. 

It is interesting to see if the interaction in eq. (33) gives a reasonable spectrum of 
baryons containing a heavy quark when the effects of operators with more than one 
derivative are neglected. The first column of energies in Table 1 gives the eigenstates 
and eigenvalues of Vj(0) in the truncated basis. Only the |0, 1/2,0), |1, 1/2, 1) and 
|1, 3/2, 1) states arc bound if g is taken positive. For the case Q = c these states have 
the right /, s, S£ and parity quantum numbers to be the Ac, Ec and S* respectively. 

■k More correctly this gives the matrix elements of the potential operator in eq. (35) divided by 
the normalization of states. The conventions for the field operators and state normalizations 
are the same as those in Ref. [14]. 
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In the large Nc limit, the N and A are degenerate and the space of states should 
be enlarged to include products of A-baryons with the P' and P* mesons. Such 
states can have the same quantum numbers as the Eg and the E*. Because the spin- 
spin interaction is suppressed by 1/iVc, including the A causes the Eg and E* to be 
degenerate with the Ac and gives the energies presented in the last column of Table 
1.-^^ The real world is intermediate between the A^"c ^ oo limit with a degenerate 
and A, and the limit where the A is neglected because it is considered to be much 
heavier than the nucleon. Including the A with a finite A — N mass difference will lead 
to a Eg state which is heavier than the Ag with a binding energy that is intermediate 
between the values in the two columns of Table 1. Thus the interaction in eq. (33) 
gives a reasonable qualitative description of heavy baryons. However it is important 
to remember that anti-baryon-meson bound states are described by taking F — > —F 
so the interaction in eq. (33) gives rise to exotic states. These states can be removed 
by including the effects of higher derivative terms in the chiral lagrangian. This is not 
neccessarily in conflict with the derivative expansion, since the normal state binding 
energy is three times the exotic state binding energy. 

In the Nc ^ oo and mg — > oo limit, the most general potential at the origin 
including all possible higher dimension operators has the form 

V = Vo{F) - lviiF){S,' - 3/2) (/2 - 3/2) , (39) 

on the truncated meson-nucleon subspace. In eq. (39) Vq and Vi are functionals of 
the soliton shape function -F(|a;|). The coefficients Vq and Vi have no definite sym- 
metry under F — > —F, so that the meson-nucleon and meson-antinucleon scattering 
potentials are not related to each other. 

The work presented in this paper is similar in spirit to that of Callan and 
Klebanov.^° The main difference is that constraints imposed by heavy quark sym- 
metry are taken into account. The bound state approach of Callan and Klebanov 
cannot be taken over to baryons containing a heavy charm or bottom quark without 
imposing heavy quark symmetry (e.g., there is no reason to expect eq. (6.2) of Ref. 
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[10] to hold). The mechanism for the parity flip presented in this paper is the same 
as that of CaUan and Klebanov. For example, if the kinetic energy of the nucleon is 
included then the ground state is an L = partial wave* instead of a Dirac delta func- 
tion. (Some of the excited negative parity baryons containing a heavy quark appear 
as bound states with odd orbital angular momentum.) However, since the unprimed 
and primed heavy meson fields are related (near the origin) by a factor of a; ■ r the 
S'-wave primed heavy meson-nucleon bound state wave functions get multiplied by 
a linear combination of Yi^m spherical harmonics in the unprimed basis and hence 
appear as P-waves. Similar changes in quantum numbers also occur in bound states 
involving monopoles.^^ 

Mark B. Wise thanks N. Seiberg for suggesting that it would be interesting to 
examine how baryons containing a heavy quark arise as solitons. 



Table 1 



States 


Energies neglecting the A 


Energies including the A 


|0, 1/2,0) 


-3^F'(0)/2 


-3gF'{0)/2 


|1, 3/2, 1),|1, 1/2,1) 


-gF'{0)/6 


-SgF'{0)/2 


|1, 1/2,0) 


gF'{0)/2 


gF'{0)/2 


|0, 3/2, 1),|0, 1/2,1) 


gF\0)/2 


gF'iO)/2 



■k With a radial wave function that is strongly peaked about = 0. 
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